Properties of Quantum Systems via Diagonalization 
of Transition Amplitudes I: Discretization Effects 
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We analyze the method for calculation of properties of non-relativistic quantum systems based on 
exact diagonalization of space-discretized short-time evolution operators. In this paper we present a 
detailed analysis of the errors associated with space discretization. Approaches using direct diago- 
nalization of real-space discretized Hamiltonians lead to polynomial errors in discretization spacing 
A. Here we show that the method based on the diagonalization of the short-time evolution oper- 
ators leads to substantially smaller discretization errors, vanishing exponentially with 1/A 2 . As a 
result, the presented calculation scheme is particularly well suited for numerical studies of few-body 
quantum systems. The analytically derived discretization errors estimates are numerically shown to 
hold for several models. In the followup paper [lj] we present and analyze substantial improvements 
that result from the merger of this approach with the recently introduced effective-action scheme 
for high-precision calculation of short-time propagation. 

PACS numbers: 02.60.-x, 03.65.-w, 31.15.X-, 71.15.Qe 

Keywords: Discretization, Exact diagonalization, Short-time propagator, Higher-order eigenstates 



I. INTRODUCTION 

In the standard operator formulation of quantum me- 
chanics, the description of a physical system is based on 
constructing the Hamiltonian operator H. Properties of 
quantum systems are obtained by solving the correspond- 
ing Schrodinger equation, 



H\ib) = E\ 



(1) 



Exact solutions can be found only for a very limited set 
of simple models. A wide variety of analytical approx- 
imation techniques has been developed in the past for 
treatment of such problems. In addition, the last two 
decades have seen a rapid growth in the application of 
different numerical methods for solving the Schrodinger 
equation. In this paper we focus on approaches based on 
real-space discretization, which usually start from some 
given finite-difference prescription. Such methods have 
been extensively studied in the past, and the main diffi- 
culties follow from the finite-difference representations of 
the kinetic operator. 

A numerical approach based on diagonalizing of the 
evolution operator, introduced in Ref. [2[, does not suf- 
fer from problems with the representation of differential 
operators on real-space grids, and has substantial ad- 
vantages in practical applications to few-body problems. 
Effectively, in this way the problem is transfered from 
that of representing the kinetic operator on a real-space 
grid to the calculating of corresponding transition ampli- 
tudes. Detailed analysis of the errors associated with the 
implementation of this approach has not been presented 
before, and is the main result of the current paper. It 
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provides full understanding of the method and allows its 
optimal use, as well as further significant improvements 
within a generalized calculation scheme, presented in our 
followup paper 

The advantages of the method discussed in this paper 
0,0, 0) HI follow from two key properties. First, the ob- 
jects being diagonalized are transition amplitudes, which 
are well defined irrespective of discretization scheme, i.e. 
the exponential of the Hamiltonian effectively regularizes 
the kinetic operator, making possible representations of 
the evolution operator that do not depend on the space 
grid. Second, the successful diagonalization of the evo- 
lution operator exp(-tH) for any time of propagation t 
immediately gives the solution of the eigenproblem for 
the Hamiltonian. Thus, the time of propagation in this 
approach is just an auxiliary parameter. Said another 
way, we use the time-dependent evolution operator to ex- 
tract time-independent information regarding the quan- 
tum system. If one could calculate transition amplitudes 
exactly, then the obtained results for the energy eigen- 
problem would not depend on the time of propagation. 
However, in practical applications one uses some approx- 
imation scheme to calculate the amplitudes, and in this 
case the precision of the obtained results for energy eigen- 
values and eigenstates does depend on time t. The gen- 
eral applicability of the outlined method follows from the 
fact that one can use short-time propagation amplitudes 
to obtain high accuracy results. 

In order to complete this numerical method and make 
it generally applicable, it is necessary to address the fol- 
lowing key questions: 

1. How to analytically estimate the effects of spatial 
discretization? 

2. How to optimize the choice of evolution time t, so 
as to minimize errors? 

3. How to accurately calculate transition amplitudes? 
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The authors in Ref. [2| have only briefly commented on 
the first two questions, and numerically determined the 
values of parameters that can be used for precise calcu- 
lations of energy eigenvalues and eigenstates for several 
models. To numerically calculate transition amplitudes, 
they exclusively relied on the naive short-time approxi- 
mation formula, 

(x\e- tH \y) « -== e"— 4 = , (2) 

V 27TI 

which is correct only to order 0(t) for energy eigenvalues. 

In this paper we address the above questions, which 
have not been fully answered before. In Section |TT] 
we present the method and notation, and identify the 
sources of the errors present in real-space discretization 
approaches. In Section [Ml we analyze in detail questions 
I and 2, and discuss the effects of discretization on the 
obtained properties of physical systems. We analytically 
derive estimates for errors coming from space discretiza- 
tion coarseness, finite size effects, and choice of evolution 
time parameter t. All the analytically derived results are 
numerically verified to hold on several instructive models. 

The followup paper [l[ continues this investigation and 
significantly improves the method by applying the re- 
cently introduced effective action approach [y, UH, d, E3] 
to completely resolve the problem formulated in question 
3. This has been addressed recently [ill E2> EH using 
various approaches. We stress that use of higher-order 
effective actions represents an efficient and numerically 
inexpensive way to calculate transition amplitudes that 
lead to many orders of magnitude increase in precision of 
calculated properties of the system. 

The expressions written throughout this paper are, 
for compactness of notation, for one particle in one di- 
mension. Extension to more particles and dimensions is 
straightforward, just as with the above short-time tran- 
sition amplitude. Note that we are working in imaginary 
time, which is well suited for numerical calculations and 
does not affect in any way energy levels nor other time- 
independent properties of the system. We have also set 
h to unity. 

II. SPACE-DISCRETIZED SCHRODINGER 
EQUATION 

In the coordinate representation the time-independent 
Schrodinger's equation takes the form 

[dy(x\H\y)(y\tl>)=E(x\rl>). (3) 

The standard way to numerically implement exact diago- 
nalization is to go from continuous coordinates x to ones 
living on a discrete space grid x n = nA, where A is a 
given spacing and n€l Integrations in the above equa- 
tion are performed using the simple rectangular quadra- 
ture rule, or some higher-order finite-difference formula. 



This completes the transition to the space-discretized 
counterpart of the continuous theory, however, to rep- 
resent this on a computer we still have to restrict the 
integers n to a finite range. This is equivalent to intro- 
ducing a space cutoff L, or putting the system in a in- 
finitely high potential box. For example, the rectangular 
quadrature rule leads to the following space-discretized 
Schrodinger equation 

N-i 

Y, H nm (mA\^) = E(A, L) (nA|^> , (4) 

m=—N 

where H nm — A ■ (nA\H\mA), N = [L/A], and square 
brackets represent the integer part of the argument. As 
a result, we have obtained a 2N x 2N matrix that repre- 
sents the Hamiltonian of the system. The eigenvalues of 
this matrix depend on the two parameters introduced in 
the above discretization process: cutoff L and discretiza- 
tion step A. Continuous physical quantities are recovered 
in the limit L — > oo and A — > 0. The outlined procedure 
is very useful in dealing with spatially localized physi- 
cal problems, such as electronic structure calculations in 
semiconductor and polymer physics (l4| . 

The two approximations involved in the discretization 
procedure, characterized by parameters A and L, are 
common steps in solving eigenproblems of Hamiltonians, 
and as such have been extensively analyzed. The im- 
posed constraint on the values of spatial coordinates to 
the finite interval (-L, L) is a valid approach for cap- 
turing information on localized eigenstates. In this ap- 
proximation the system is effectively surrounded by an 
infinitely high wall, and as the cutoff L tends to infinity, 
we approach the exact energy levels always from above 
[TEl llq ], which is a typical variational behavior. There- 
fore, we designate errors associated with the cutoff L as 
variational. The effects of the discretization step A are 
much more complex, and follow from the fact that the 
kinetic energy operator cannot be exactly represented on 
finite real-space grids. For example, a typical naive dis- 
cretization of the kinetic energy operator gives in our 
notation the following Hamiltonian matrix elements [l7| 

( l/A 2 + y(?iA) ifn = m 
H nm = \ -1/(2A 2 ) if|n-m| = l (5) 
I otherwise. 

Note that in the absence of a potential term V in the 
Hamiltonian, the above definition corresponds to a tight- 
binding model [l7j . This prescription leads to numerical 
results for eigenvalues which in the A — > limit converge 
to the exact continuum values as A 2 . The errors asso- 
ciated with this approach have non- variational behavior, 
i.e. the obtained results are not always upper bounds of 
the exact energy levels. Several papers discuss this issue 
and analyze the behavior of errors in the direct diagonal- 
ization approach (for more details, see Refs. [HI, E|| [2(| 
and references therein). The state-of-the-art in this ap- 
proach is a set of systematically improved prescriptions 
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FIG. 1: Harmonic oscillator transition amplitude as a func- 
tion of coordinates x and y for t = 1, frequency u) = 1 and 
mass M = 1. 

for discretization of the kinetic energy operator, which 
speeds up convergence to the continuum limit to higher 
powers of A 2 . However, within this approach conver- 
gence is always polynomial in A. Some recent results 
[19L [20l | also exist on extensions of this approach that 
provide effective variational behavior of the discretized 
kinetic energy operator. 

As outlined in the Introduction, in this paper we focus 
on an alternative approach, based on solving the eigen- 
problem of the corresponding transition amplitudes as 
proposed in The central equation is 

AT-l 

]T A nm (t) (mA\i>) = e- tE ^> L ^ <nA|^> , (6) 

ro=— N 

where A nm (t) = A ■ A(nA, mA; t) = A • (nA|e _t& |mA). 
In this approach the time of evolution t plays the role 
of an auxiliary parameter. This parameter is not re- 
lated to the discretization, and in a continuous theory it 
does not affect the obtained eigenvalues and eigenstates. 
However, in a discretized theory the numerically calcu- 
lated eigenvalues and eigenstates will necessarily depend 
on this parameter as well, as emphasized by the right- 
hand size of Eq. Therefore, the original problem 
is now transformed into the eigenproblem of the matrix 
A nm (t) , whose indices take all integer values in the range 
-N <n,m<N, where N = [L/A]. 

Fig.[T]shows how a typical transition amplitude, in this 
case that of a harmonic oscillator, depends on coordi- 
nates x and y. As can be seen from the figure, transition 
amplitudes are spatially well localized. This is particu- 
larly simple to understand for the short times of prop- 
agation that we consider. In this case the kinetic term 
exponentially localizes the transition amplitude matrix 
to the vicinity of the main diagonal. Similarly, the po- 
tential brings about exponential localization along the 
main diagonal around its minimum. The localization of 
dominant values of the transition amplitude to a small 
area in the x — y plane gives practical justification for 
introduction of space cutoff L in this approach. 



In continuum theory, the transition amplitude eigen- 
problem is mathematically equivalent to the Schrodinger 
equation. It is important to stress, however, that the 
procedure of space discretization introduces important 
differences between eigenproblems Q and (|6]). In par- 
ticular, as we will show in the next section, the proce- 
dure based on the diagonalization of transition ampli- 
tudes leads to much faster (non-polynomial) convergence. 
An illustration of the relation of these two calculation 
schemes is shown in Fig. [5] which compares the exact 
parabolic dispersion of a free particle in a box with nu- 
merical calculations based on the diagonalizations of the 
Hamiltonian and of the transition amplitudes. From the 
figure we see that the time parameter t in the transition 
amplitude approach plays an important role. Increase of 
t gives better agreement with the exact dispersion rela- 
tion. 
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FIG. 2: (Color online) Eigenspectrum of a free particle in a 
box. Eigenvalues Ek are given as a function of level num- 
ber k. The solid line gives the exact parabolic dispersion 
Ek = 7r 2 (fc + 1) 2 /8L 2 , while the dashed line presents results 
calculated in the tight-binding approximation. The graph also 
shows numerical results obtained by the diagonalization of 
transition amplitudes for different values of time of evolution 
t. All the numerical calculations are for L — 6 and A = 0.25, 
hence N = L/A = 24. 



III. DIAGONALIZATION OF 
SPACE-DISCRETIZED TRANSITION 
AMPLITUDES 

The free particle transition amplitude 

^■« t >-v53-'(- fc ir ! ) (7) 

satisfies 

J dxAf ree (x,y 7 t) = l . (8) 
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The consequence of this is conservation of probability. 
In the space-discretized analogue of this model x = nA, 
y = to A, and the transition amplitude is A/jj^ e (i) = 
A Af ree (nA, toA; t). Using the Poisson summation for- 
mula 



E ex p (- an2 ) = E ex p (~ 



(9) 



we find that the space discretized free particle amplitude 
satisfies 



E A free 
nm 



(*) 



, 2tt 2 
l + 2exp ( —gft 



Conservation of probability is thus obtained only in the 
continuum limit A — > 0. Note that the effect of dis- 
cretization is non-perturbative in discretization step A, 
i.e. it is smaller than any power of A. The effect of dis- 
cretization is also universal in that it holds for all models, 
since the free particle transition amplitude is the domi- 
nant term in the short time expansion of the transition 
amplitude of a general theory. 

To show this explicitly we use the short time expansion 
of the transition amplitude of a general theory [6| to show 
that 

J dxA(x,y : t) = — 1= J dxe-^ x2 ^2t l fi(x,y) , 

(11) 

where /o = 1 and the other /; are given functions of the 
potential and its derivatives. Writing the even part of 
fl(x,y) as g(x 2 ,y) we find 



da; A(x, y; t) 



(12) 



Similarly, using the above Poisson summation formula, 
we find 



^ A nm (t) - / dxA(x,y;t) 



(13) 



This universal and non-perturbatively small deviation 
from the continuum indicates that one should center 
numerical calculation schemes on transition amplitudes 
rather than the Hamiltonian. By diagonalizing the tran- 
sition amplitude for any time of propagation t we obtain 
the energy eigenvalues and eigenfunctions 



dy A(x,y;t)ip k (y) 



-tE k 



(15) 



To solve this numerically we first discretize space with 
discretization step A, and second we introduce a spatial 
cut-off L such that |a;| < L. Amplitudes are now 2Nx2N 
matrices whose diagonalization leads to 2N eigenstates 



-J 




Performing the indicated differentiations the right hand 
side becomes exp(— ^-t) • Y^i hl{v)t l ■ O ne could now 
calculate the hi from the short time expansions /;. The 
p-level effective action gives a short time expansion that 
is truncated at order t p . As a result J2i hi{y)t l is a poly- 
nomial in time of order p. The dominant short time be- 
havior is thus given by the universal exponential term. 
As a result the transition of a general model to its space 
discretized form is given by 



^ A nm (t) - / dx A(x, y; t) ~ exp 



2^2 
' A 2 



(14) 



FIG. 3: (Color online) (a) Plot of \E (A,L,t) - E \ for a 
free particle in a box as a function of A for different values 
of time of evolution t and L = 6. For comparison, we also 
plot the corresponding deviations of numerical results (des- 
ignated by H) obtained using direct diagonalization of the 
space-discretized Hamiltonian, denned by Eq. ([5}. The inset 
gives a closer view of direct Hamiltonian diagonalization er- 
rors, since they have much weaker dependence on the spacing 
A. (b) This plot shows how the deviations |£^(A, L, t) — Ej,\ 
depend on t for several energy levels k. The parameters used 
are L = 6, A = 0.2. In both plots the dashed lines represent 
discretization error estimates given in Eq. (|17|l . 
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ipk and eigenvalues e -* E k(A,L,t) _ 

As we have seen, discretization introduces a non- 
perturbatively small error in transition amplitudes pro- 
portional to exp (— 27r 2 £/A 2 ). We should therefore ex- 
pect the discretization error for energy eigenvalues to be 

£ fe (A,L,i)-£ fe ~-iexp(-^ r ^ . (16) 

We have numerically investigated this for a diverse set of 
models and have shown the above relation to hold in all 
cases. It is also illustrative to verify this for analytically 
tractable models. Using the known analytical expres- 
sions for transition amplitude and energy eigenstates for 
a free particle in a box [2l|, [22j, as well as the Poisson 
summation formula in Eq. ©, we find that the energy 
eigenstates of the space discretized model satisfy 

2 2„2 /7T 2 (fc+l)A , , 

£ fc (A,L,t)-£ fc = --e-^ f cosh( \ A ' J , (17) 

where Ek — - and k = 0, 1, 2, . . . As expected, 

the universal term gives the dominant A dependancc. 
One obtains similar analytical results for the case of the 
harmonic oscillator. 

The non-perturbatively small effect of spatial dis- 
cretization is the reason why the new method highly out- 
performs direct diagonalization of the Hamiltonian and 
leads to much smaller errors for the same size of dis- 
cretization step A. In addition to this the free param- 
eter associated with the method, the time of evolution 
t, can be used to further minimize errors. As illustrated 
in Fig. [2l while keeping A fixed, we can adjust time t 
to obtain much smaller errors and practically reproduce 
the exact spectrum of the theory. This is also evident in 
Figs. [3^ and[3]3, where we see that by adjusting t, errors 
can be reduced by orders of magnitude for fixed value of 
discretization step A. 

We next consider a harmonic oscillator. Fig Q] shows 
how the presented method may be used to obtain en- 
ergy eigenvalues to high levels. The numerical calcula- 
tions agree well with the well known linear dispersion 
of the harmonic oscillator. Figs. [5^, and b display re- 
spectively the A and t dependance of the deviations 
\Ek(A, L, t) — Ek\, showing agreement with the analyti- 
cally derived estimate of the discretization error given in 
Eq. p6[) . In order to achieve such a high accuracy of nu- 
merical results as presented on all grap hs, we have used 
the Mathematica software package [231 ]. 

Fig. 03 shows how the deviations \E k (A, L,t) — Ef.\ 
depend on t for several levels k. The plot corresponds to 
the harmonic oscillator but is typical of a general theory. 
The saturation of errors for large t comes about when 
the discretization error, given by the universal estimate 
in Eq. p^|) . becomes smaller than the error due to space 
cutoff L. Analytical estimates for cutoff error are given 
at the end of this section. At this point we just mention 
that the finite size effects can already be seen in Fig. 0] 
where for high values of level number k numerical results 



start to move away from the linear dispersion character- 
istic of a harmonic oscillator to the parabolic dispersion 
characteristic of a box potential. 

So far we have considered only integrablc models, i.e. 
models for which we know the exact transition ampli- 
tudes. As a result we have thus far encountered and an- 
alyzed only two sources of errors: those associated with 
discretization step A and cutoff L. The vast majority of 
models are not integrable. The outlined method is still 
applicable if one uses some approximation for calculat- 
ing transition amplitudes. In a previous series of papers 
we have used the method of effective actions to calcu- 
late short time expansions of transition amplitudes of a 
general theory to high order p. For the case of a general 
many particle theory in arbitrary number of dimensions 
we have obtained closed form expressions for expansions 
up to p — 10. The analytical procedure is substantially 
simplified for certain potentials. In particular, for poly- 
nomial interactions we have obtained expressions to level 
p = 144. This high level of precision makes these short 
time expansion formulas ideal for use in the method out- 
lined in this paper. Still, the diagonalization of approx- 
imate transition amplitudes introduces a third source of 
error proportional to t p . 

The overall error is minimized when all three sources 
of error are approximately equal. The key point is that 
we have simple analytical estimates for all three errors. 
The universal behavior of the discretization error sub- 
stantially simplifies the process by which one chooses the 
values of parameters A, L, and t that minimize the over- 
all error. The details of this are presented in the followup 
paper. 

Fig. [5] displays \Eq(A, L,t) — Eq\ as a function of dis- 
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FIG. 4: (Color online) Harmonic oscillator dispersion relation. 
The solid line gives the exact linear dispersion Ek = k + 
1/2. The points correspond to numerically calculated energy 
eigenvalues Ek as function of level k. We show the results of 
the diagonalization of transition amplitudes for several values 
of t. In this plot L = 12, A = 0.25, the frequency is u = 1, 
mass M = 1. 
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FIG. 6: (Color online) Plot of \E Q [A, L, t) - E \ for an anhar- 
monic oscillator with potential V = \Mlo 2 x 2 + ^x 4 given as 
a function of A for different values of time of evolution t. The 
parameters used are L — 6, harmonic frequency uj = 1, mass 
M — 1, and anharmonicity g — 48. Transition amplitude ma- 
trix elements were calculated using p = 18 effective actions 
[(J . The high precision value for the exact ground energy that 
we compare to was calculated in Ref. [24| . Dashed lines cor- 
respond to the discretization error in Eq. (jXTJ) . For compar- 
ison, we also plot the corresponding deviations of numerical 
results (designated by H) obtained using direct diagonaliza- 
tion of the corresponding space-discretized Hamiltonian. The 
inset gives a closer view of direct Hamiltonian diagonaliza- 
tion errors, since they have much weaker dependence on the 
discretization step A. 



FIG. 5: (Color online) (a) Plot of \E (A,L,t) - E \ for a 
harmonic oscillator as a function of discretization step A for 
different values of time of evolution t with L — 12, u — 1, and 
M — 1. For a comparison, we also plot the corresponding re- 
sults (designated by H) obtained using direct diagonalization 
of the space-discretized harmonic oscillator Hamiltonian. The 
inset gives a closer view of direct Hamiltonian results, since 
they have much weaker dependence on the discretization step 
A. (b) Plot of the deviations \Ek(A, L,t) — Ek\ given as a 
function of time t for several levels k. The parameters used 
are L — 12, A = 0.1, u) = 1, M = 1. The observed saturation 
of errors for large t comes about when the discretization error 
becomes smaller than the error due to space cutoff L. In both 
plots the dashed lines correspond to the discretization error 
estimate for Eq given in Eq. (|17[) . 



cretization step A for the case of an anharmonic oscil- 
lator with potential V — ^Mlo 2 x 2 + ^x 4 . The param- 
eters used in the plot are L — 6, u> — 1, M = 1, and 
anharmonicity g — 48. The transition amplitude ma- 
trix elements were calculated using p — 18 effective ac- 
tions @. The high precision value for the exact ground 
energy that we compare to was calculated in Ref. [24| . 
As we can see, even though we are dealing with a rela- 
tively strong anharmonicity, the numerically calculated 
values stay right on the dashed lines corresponding to 



the universal discretization error just as in the case of 
the previously considered integrable models. This is in 
complete agreement with our analytical derivation of the 
discretization error. In the followup paper we numeri- 
cally investigate a variety of different interacting models 
and in all cases document the validity of this formula. 

As can be seen from Fig. [SJ the numerical results 
clearly demonstrate that the A-dependence of errors 
within our calculation scheme highly outperforms the 
polynomial dependence in A 2 obtained by the direct di- 
agonalization of the Hamiltonian. This is true even for 
short times of propagation t. Although interaction terms 
in the potential affect the numerical values of errors, di- 
agonalization of the transition amplitudes still substan- 
tially outperforms diagonalization of the Hamiltonian, 
and is the preferred method. This success is a conse- 
quence of the non-perturbative behavior of the spatial 
disretization error within this calculation scheme. This 
leads us to the key conclusion that discretization param- 
eters can be always optimized so that presented approach 
of solving eigenvalue problem of space-discretized tran- 
sition amplitudes highly outperforms direct diagonaliza- 
tion of the space-discretized Hamiltonian. The contin- 
uum limit A — > is far more easily approached in the 
first case and the corresponding discretization errors are 
substantially smaller for the same discretization coarse- 
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ness. From the numerical point of view, as the value of 
parameter A directly determines the size of the matrix 
to be diagonalized, the computational cost for the same 
precision is significantly reduced. 

We end by looking at finite size effects, i.e. errors re- 
lated to introduction of space cutoff L. For any theory 
with non-trivial potential, the cutoff L is artificially in- 
troduced and it affects the obtained energy eigenvalues, 
as we have already discussed. To estimate the effects 
of the cutoff, we first note that they are closely related 
to the spatial extent of the potential V, as well as the 
spatial extent of eigenfunctions of the system: errors in 
the corresponding energy eigenvalues can be considered 
small only if the eigenstates ipk(x) are well localized in 
the interval |x| < L. 

The effects of space cutoffs have beenpreviously stud- 
ied for continuous-space theories [HI, [l6j]. The shift in 
energy level Ek(L) — Ek is found to be positive in this 
case, and approximately given by the formula 

E '"'- ,i - aw (i ») ' (18) 

where a is an appropriately chosen value of coordinate x 
such that it is larger than and well away from the largest 
zero of ipk(x) but smaller than and well away from the 
space cutoff L. The constant Cf-(a) depends on the nor- 
malization of eigenfunction and the choice of parameter 
a. For example, the ground state has no zeros, and we 
can always choose the value a = 0. In that case, constant 
C o (0) is given by 

Co(0)= (^J L dx\Mx)\^J , (19) 

where we assume that the eigenfunction tpo(%) is normal- 
ized as usual, dx \ipo(x)\ 2 = 1. 

In practical applications, when we use diagonalization 
of the discretized transition amplitudes, the errors in en- 
ergy level will necessarily also depend on the parameter t 
and other discretization parameters. Here we give a sim- 
ple estimate of ground energy errors that follows from 
the spectral decomposition of diagonal amplitudes. For 
large t we have A(x,x;t) ~ \^po{x)\ 2 e~ Eot . Integrating 
this we find an approximate result for the ground energy 
of a system with cutoff L 

E (L,t) --\n [ dxA(x,x;t), (20) 

t J-L 

In the L — > oo limit we recover the exact ground energy, 
so that a simple estimate of finite size effects on E is 
given by 

E Q {L,t)-E Q &\ f dx\Mx)\ 2 - (21) 

1 J\x\>L 

Although the above equation is just a rough estimate of 
the errors introduced by a space cutoff L, Fig. [7] shows 
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FIG. 7: (Color online) Deviations £*(A, L, t) — Ek for a har- 
monic oscillator as a function of space cutoff L for differ- 
ent values of time of evolution t. The parameters used are 
A = 0.1, w = 1, M = 1. Solid thin lines give the dominant 
behavior of Eq. (|21l) . The dashed thick lines correspond to 
the error estimate in Eq. (|18|) . 

that it is in good agreement with numerical results for 
the harmonic oscillator. In order to clearly demonstrate 
L-dependence of errors in this graph, we have used small 
value of the discretization step A, such that discretization 
errors can be neglected. The dashed line in the figure 
represents error estimates given by Eq. ()18|) . 

Using the data from Fig.[7]we can now fully explain the 
saturation of errors observed in Fig. O). The value of the 
cutoff L used to obtain this data was L — 12. As can be 
seen from Fig. [3 this value of the cutoff parameter yields 
an error of the order 10~ 65 for the ground-state energy 
for t ~ 0.1, and of the order 10 -40 for energy eigenlevel 
En. These values exactly correspond to the saturated 
errors in Fig. [5Jd. 

Although in the general case the eigenstates that come 
into Eqs. (fT5| and (|2"Tj) are not known, we can still use 
them in conjunction with other approximation techniques 
to estimate finite size effects. We also see that, due to 
the larger spatial extent of higher energy eigenstates, the 
cutoff-related errors are minimal for the ground energy. 
Note however that one is not really interested in the pre- 
cise calculation of finite size errors, but only needs to 
estimate the minimal size of the cutoff L for which finite 
size effects are negligible. For that purpose one can use 
either of the above approximate formulas. 



IV. CONCLUSIONS 

The current paper is the first in a scries of two publi- 
cations dealing with the properties of quantum systems 
calculated from the the diagonalization of transition am- 
plitudes. In this paper we have focused on analyzing the 
errors associated with real-space discretization and finite 
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size effects. In particular, we have shown that within this 
calculation scheme spatial discretization leads to a uni- 
versal and non-perturbatively small discretization error. 
This highly outperforms the usual polynomial behavior 
of errors in approaches based on the diagonalization of 
space-discretized Hamiltonians. In addition to providing 
a full understanding of the numerical method based on 
diagonalization of the evolution operator, we have also 
derived analytical estimates for all the errors involved 
within this approach. In practical applications, the de- 
rived analytical results make it possible to optimize pa- 
rameters of the method so as to minimize errors in cal- 
culated energy eigenvalues and eigenstates in the case of 
a general theory. 

The second paper in the series [1] builds on these re- 
sults, extending them through systematic improvement 
of short-time propagation using the effective action ap- 
proach [6( ■ This effectively solves the problem of the ac- 



curate calculation of evolution operator matrix elements 
and significantly reduces errors related to the time of evo- 
lution parameter. 
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